arXiv:1502.01750v3 [math.PR] 12 Feb 2015 


October 17, 2014; To appear in: Advances in Applied Probability 


GAUSSIAN RANDOM PARTICLES WITH FLEXIBLE HAUSDORFF DIMENSION 

LINDA V. HANSEN,* Varde College 

I 

THORDIS L. THORARINSDOTTIR,** Norwegian Computing Center 
EVGENI OVCHAROV,*** Heidelberg Institute for Theoretical Studies 
! TILMANN GNEITING,**** Heidelberg InsUtute for TheoreUcal Stupes and Karlsruhe InsUtute of Technology 
^ DONALD RICHARDS,***** Penn State University 

Abstract 

Gaussian particles provide a flexible framework for modelling and simulating 
three-dimensional star-shaped random sets. In our framework, the radial 
function of the particle arises from a kernel smoothing, and is associated with 
an isotropic random field on the sphere. If the kernel is a von Mises-Fisher 
density, or uniform on a spherical cap, the correlation function of the associated 
random field admits a closed form expression. The Hausdorff dimension of the 
surface of the Gaussian particle reflects the decay of the correlation function at 
the origin, as quantified by the fractal index. Under power kernels we obtain 
particles with boundaries of any Hausdorff dimension between 2 and 3. 
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1. Introduction 


Mathematical models for three-dimensional particles have received great interest in astronomy, botany, geology, 
material science, and zoology, among many other disciplines. While some particles such as cryst als have a rig id 


shape, many real-world objects are star-shaped, highly structured, and stochastically varying (IWicksell . 


1925 


Stovan and Stovan 


in high demand. 


1994^. As a result, flexible yet parsimonious models for star-shaped random sets have been 


Gren ander and Miller (119941) proposed a model for two-dimensional featureless objects with no 


obvious landmarks, which are represented by a deformed polygon along with a Gaussian shape model. This was 
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inv estigated furth e r by 


by 


Hobolth et al 


Kent et al 


(120031) . 


Miller et al 


harmonics and was studied by 


(200C) and 


Hobolth et al 


(12002h . and a non-Gaussian extension was suggested 


19941 ) proposed an isotropic deformation model that relies on spherical 


Hobolth (2003 ), where i t was applied to m onitor tumour growt h. A related Gaussian 


random shape model was studied by 
desert dust particles. 


Muinonen et al 


(1996|) and used by 


Munoz et al. 


(|2007l) to represent Saharan 


In this paper we propose a flexible framework for modelling three-dimensional star-shaped particles, where the 
radial function is a random field on the sphere that arises through a kernel smoothing. Specifically, let Y C R 3 be a 
three-dimensional compact set, which is star-shaped with respect to an interior point o. Then there is a one-to-one 
correspondence between the set Y and its radial function X = (A(u) : u £ S 2 }, where 


X(u) = max{r > 0 : o + ru £ Y}, u £ S 2 , 


with S 2 = (i £ I 3 : ||a;|| = 1} denoting the unit sphere in R 3 . We model X as a real-valued random field on § 2 via 
a kernel smoothing of a Gaussian measure, in that 

X(u) = [ K(v,u)L(dv), u£E> 2 , (1) 

Js 2 

where K : 8 2 x § 2 —> R is a suitable kernel function, and L is a Gaussian measure on the Borel subsets of § 2 . That 
is, L(A) ~ A/”(/i A(A), cr 2 A(A)) with parameters /j £ R and cr 2 > 0, where A(A) denotes the surface measure of a 
Borel set A C § 2 , with A(§ 2 ) = 47 t. 

If X were a nonnegative process, the random particle could be described as the set 


Y = |^J {o + ru : 0 < r < X(u)} C R 3 , 

u £§ 2 

so that the particle contains the centre o, which without loss of generality can be assumed to be the origin, and the 
distance in direction u from o to the particle boundary is given by X{u). A potentially modified particle Y c arises 
in the case of a general, not necessarily nonnegative process, where we replace X(u) by X c (u) = max(c, X(u)) for 
some c > 0. We call Y or Y c a Gaussian particle, with realisations being illustrated in Fi gure [H The Ga ussian 


particle framework is a special case of the linear spatio-temporal Levy model proposed by 


Jonsdottir et al. 


( 20081 ) 


in the context of tumour g rowth. A lterna t ively, it can be seen as a generalisation and a three-dimensional extension 
of the model pro posed by 


models of 


Miller et al 


Hobolth et al 


1 19941) and 


(2003^, w hile a lso being a generalisation of the Gaussian random shape 


Muinonen et al 


(1996)). 


The realisations in Figure |T] demonstrate that the boundary or surface of a Gaussian particle allows for regular 
as well as irregular behaviour. The roughness or smoothness of the surface in the limit as the observational scale 
becomes infinitesimally fine can be quantified by the Hausdorff dimension, which for a surface in R 3 varies between 
2 and 3, with the lower limit corresponding to a sm ooth, differentiable surface, and the upper limit corresponding to 


an excessively rough, space-filling surface ((Falconer, 


1990I ). The Hausdorff dimension of the surface of an isotropic 


Gaussian particle is determined solely by the behaviour of the correlation function of the associated random field 
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Figure 1: Gaussian particles with mean px = 100 and variance ax = 10, using the power kernel Q with q = 0.05 (left), 
q = 0.25 (middle) and q = 0.5 (right). The Hausdorff dimension of the particle surface equals 2 + q. 


on the sphere. We investigate the properties of Gaussian particles under parametric families of is otrop ic kernel 


functions, including power kernels, and kernels that are proportional to von Mises-Fisher densities ( Fishe r et al 


19871) . or uniform on spherical caps. Under power kernels we obtain particles with boundaries of any Hausdorff 


dimension between 2 and 3. Von Mises-Fisher and uniform kernels generate Gaussian particles with boundaries of 
Hausdorff dimension 2 and 2.5, respectively. 


The remainder of the paper is organised as follows. Section [2] recalls basic properties of the radial function in the 
Gaussian particle model ©■ In Section [3] we show how to derive the Hausdorff dimension of an isotropic Gaussian 
particle from the infinitesimal behaviour of the correlation function of the underlying random field at the origin. 
Section [I] introduces the aforementioned families of isotropic kernels and discusses the properties of the associated 
correlation functions and Gaussian particles, with some technical arguments referred to an appendix. Section [5] 
presents a simulation algorithm and simulation examples, including a case study on celestial bodies and a discussion 
of planar particles. The paper ends with a discussion in Section [6] 


2. Preliminaries 

The properties of the random function m that characterises a Gaussian particle process depend on the kernel 
function K. We assume that K is isotropic, in that K(v,u) = k(d(v,u)) depends on the points v,u £ S 2 through 
their great circle distance d(v,u) £ [0,7r] only. As d(v,u) = arccos(u • i>), th is is eq uivalent to assuming that the 


Jonsdottir et al 


( 2008 ) in concert with the rotation 


kernel depends on the inner product u ■ v only. Results of 
invariance property following from an isotropic kernel imply that the mean function E(X(u)) and the variance 
function Var (X(u)) are constant, that is, 

fix = E(X(u)) = tic\ and a\ = Var(X(u)) = cr 2 ci 


for u £ § 2 , where we assume that 


c n = [ k(d(v,u)) n du 

Js 2 


is finite for n = 1,2. 
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Note that X is a stochastic process on the sphere (Jones, 1963T ). whose covariance function is given by 


Cov(A(«i), X (u 2 )) = a 2 / k{d(v,ui)) k(d{v,u 2 )) dv, «i ,«2 6 S 2 , 

J s 2 

Under an isotropic kernel, the random field X is isotropic as well, and it is readily seen that Corr(X(iii), X{u 2 )) = 
C(d(ui,u 2 )), where 


2 

C(9 ) = — / k(j]) fc(arccos(sin 9 sin tj cos (j> + cos 9 cos77)) d</> sinry d?y, 0 < 6 

C2 Jo Jo 


< 7T, 


( 2 ) 


is the correlation function of the random field X. As recently shown bv IZiegell ([2014) ), any continuous isotropic 
correlation function on a sphere admits a representation of this form. 


3. Hausdorff dimension 


uni . 


For e > 0, an e-cover of Z is a 


The Hausdorff dimension of a set Z C M“ is defined as follows (jHausdorfl . 
countable collection {Bi : i = 1, 2,...} of balls Bi C R d of diameter \Bi \ less than or equal to e that covers Z. With 


H S (Z) = lim inf { ^2\Bi\ s : {Bi : i = 1,2,...} is an e-cover of Z{ 

denoting the S -dimensional Hausdorff measure of Z, there exists a unique nonnegative number (5o such that H s (Z) = 
00 if S < <5o and H S (Z) = 0 if S > Jo- This number 5q is the Hausdorff dimension of the set Z. Note that we have 
define d the Haus dorff measure using c overings with balls. This approach is consistent with the treatments given by 
Adler (20101 and Hall and Roy (1994) and simplifies the presentation. 

As X is a kernel smoothing of a Gaussian measure, A' has Gaussian finite dimensional distributions and thus 
is a Gaussian process. While there is a wealth of results on the Hausdorff dimension of the graphs of stationary 
Gaussian random fields on Euclidean spaces, which is determine d by t he infin itesim al behaviour of the correlation 


function at the origin, as formalised by the fractal index (IHall and Rov 


1994: 


Adler 


2010f) . we are unaware of any 


extant results for the graphs of random fields on spheres, or for the surfaces of star-shaped random particles. 

We now state and prove such a result. Toward this end, we say that an isotropic random field X on the sphere 
with correlation function C : [0, tt\ —> K has fractal index a > 0 if there exists a constant b > 0 such that 


lim C(0) ~ = b. 


04.0 


9 ° 


( 3 ) 


The fractal index exists for essentially all correlation functions of practical interest, and it is always true that 
a € (0,2]. To see this, suppose that C : [0,7r] —> R is an isotropic correlation function on the two-dimensional 
sphere. Clearly, C also is an isotropic correlation function on the circle, and its even, 2-7 t periodic continuation to 
R is a station ary correlation function on the real line. Therefore, the corresponding restriction on Euclidean spaces 
(Adler, 2010 . p. 200) applies, in that a £ (0,2]. 

The following theorem relates the Hausdorff dimension of the graph of an isotropic Gaussian random field X on 
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the sphere § 2 to its fractal index. The proof employs stereographic projections that allow us to draw on classical 
results in the Euclidean case. 

Theorem 1. Let X be an isotropic Gaussian random field on § 2 with fractal index a £ (0, 2], Consider the random 
surface 

Z c = {(u,X c (u)) : u £ § 2 }, 

where X c (u ) = max(c, X(u)) with c > 0. Then with probability one either of the following alternatives holds: 

(a) If max ug §2 X(u ) < c, the realisation of Z c is the sphere with radius c and so its Hausdorff dimension is 2. 

(b) If max ug §2 X(u) > c, the realisation of Z c has Hausdorff dimension 3 — 

Proof. The claim in alternative (a) is trivial. To prove the statement in alternative (b), we assume without loss 
of gener ality that X(uo) > c, where uq = (0,0,1). The sample paths of X are continuous almost surely according 


to 


Gan golli (119671 Theorem 7.2). Thus, there exists an e £ (0, i) such that X(u) > c for u in the spherical cap 


§ 2 = {u £ § 2 : d(u,ug) < e} of radius e centred at uq. Let II : § 2 —> B c denote a stereographic projection 
that maps (0,0,1) to (0,0), where B e = {x = (xi/x^) £ K 2 : x\ + x\ < e 2 }. A stereographic projection is a 
local diffeo morphism, II thus is differentiable and has a differentiable inverse II -1 , which is locally bi-Lipscliitz 


(do Carrno 


A > 1 with 


1976). We may therefore assume that e is small enough so that for all x. x' £ B e there exists a constant 
1 


A 


II* — x'\\ < ||n -1 (a;) — II _1 (a; , )|| < A\\x — x'\\, 


(4) 


where || • || denotes the Euclidean norm on K 2 or K 3 , respectively. Without loss of generality, we may in the following 
consider conditional probabilities which depend < 

C R 2 be given by W(x) = J 5 f(II” 1 (a;)). From 


Xue and Xiao 

(2011 

. Theorem 5.1), see also Chapter 8 in 

Adler 


(120101 ). the graph GrW = {(a;, W(x)) : x £ B e } has Hausdorff dimension 3 — f almost surely if there exists a 
constant Mg > 1 such that 


jr E \*i - X T ^ - w/ ^')) 2 < Mo E l*i - 

0 j =1 3=1 

for all x,x' £ B e . Letting fi{x,x') = d(n _1 (x), n” 1 (^ / )), we have 

E(W(x) - W(x')) 2 = 2a% [(7(0) - C(&{x, *'))], 


(5) 


( 6 ) 


where (7 : [0, tt\ —> M is the correlation function of the isotropic random field X. As chord length and great circle 
distance are bi-Lipschitz equivalent metrics, there exists a constant B > 1 such that 

\ lin-H*) - n-V)|| < fi(x,x’) < B nn- 1 ^) - n~V)ll- 


(7) 
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As the random field X is of fractal index a, there exists a constant M\ > 0 such that 
2 

Xj-xX < 2 1- f \\x-x'\\ a < 2 1 ~^A a B a 9(x,x') a < [(7(0) - 

1=i 

for x, x' £ B e and e > 0 sufficiently small, where the first estimate is justified by Jensen’s inequality and the second 
by (p£|) and (|7|) . Similarly, there exists a constant A ^2 0 such that 

2 

M 2 [(7(0) - C(0(z,*'))] < - xT 

3 =1 

for all x, x' £ B e and e > 0 sufficiently small. In view of equation ©, this proves the existence of a constant Mq > 1 
such that © holds, given that e > 0 is sufficiently small. 

Now, consider the mapping ( from B ( x I to x 1 defined by ((x,r) = (II _1 (a;),r), so that £(GrVb) = 
{(u, X(u)) : u £ 8g}. The identity 


||C0r, r) C(x', r')|| 2 = \\U~\x) - If V)f + |r - r' 


J 12 


along with ([4]) implies C, to be bi-Lipschitz. Therefore by Proposition 3.3 of Falconer (I19901) . the partial surface 
(it, X(u) ) : u £ §^} has Hausdorff dimension 3 — almost surely. Invoking the countable stability property 


(Falconer 


1990L p. 49), we see that the full surface Z c = {(u,X c (u)) : u £ 8^} also has Hausdorff dimension 3 — 


almost surely. 


4. Isotropic kernels 

It is often desirable that the surface of the partic l e process possesses the same Hausdorff dimensio n as that of 


the real-world particles to be emulated (IMandelbrot 


1982 


Orford and Whallev 


1983; 


Turcotte 


198711 . With this 


in mind, we introduce and study three one-parameter families of isotropic kernels for the Gaussian particle process 
CD- The families yield interesting correlation structures, and we study the asymptotic behaviour at zero, which 
determines the Hausdorff dimension of the Gaussian particle surface. 

4.1. Von Mises—Fisher kernel 

Here, we consider k to be the unnormalised von Mises-Fisher density, 


k{9) = e a 


0 < d < 7T, 


with parameter a > 0. The von Mises-Fisher density with parameter a > 0 is widely used in the analysis of spherical 


data ( Fisher et al 


19871) . and in this context a is called the precision. Straightforward calculations show that 
2 sinh (a -\/2(T+~cosl?) 


C{9) = 


sinh(2a) ^2(1 + cos0) 


0 < 9 < 7T, 
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from which it is readily seen that the fractal index is a = 2. The surfaces of the corresponding Gaussian particles 
are smooth and have Hausdorff dimension 2, independently of the value of the parameter a £ R. 

4.2. Uniform kernel 

We now let the kernel k be uniform, in that 


k(8) = 1 {6 < r), 0 < 8 < ir, 


with cut-off parameter r £ ( 0 , ^]. 
correlation function is 


As shown in the appendix of Tovchi g recli ko and Vakser ( 20011) . the associated 


C(8) = 


1 


7 r — arccos 


cos 8 — cos 2 r 


7 r (1 — cosr) 

— 2 cos r arccos | cot r 


1 — cos 2 r 
1 — cos 6 
sin 0 


1(8 < 2r), 


0 < 8 < ir. 


In particular, if r = ^ then C(8) = 1 — — decays linearly throughout. Taylor expansions imply that the correlation 
function has fractal index a = 1 for all r € (0, |), so that the corresponding Gaussian particles have non-smooth 
boundaries of Hausdorff dimension |. 


4.3. Power kernel 

Our third example is the power kernel where the isotropic kernel k is defined as 

k(8) = - 1 , 0 < 8 < 7T, 

with power parameter q £ (0,1). The associated correlation function (121) takes the form 


where 


C(8) 


2 r 


C2 Jo 


(• 7 t 9 A 9 — l) sin A 



1) d <t> dA, 


t(8, A, <j>) = sin 8 sin A cos 0 + cos 8 cos A, a(9, A, <j>) = arccos t(8, A, 0), 


and 

A(A) = {(p € [0, 7 r] : 0 < a(8, A, 4>) < 7 r}. 
The normalising constant C 2 is here given by 


C2 = 27T 


/»7T 

/ U 9 A- 

^0 


— l) 2 sinA dA 


(g)3 y- (—l) J 7r 2j + 3 1 

6 ^ (2j + 1 )! (1 -q + j) 3 ' 


( 8 ) 


(9) 


where ( 0)3 = a(a + l)(a + 2). This expression for C 2 is obtained by expanding sin A in a Maclaurin series and then 
integrating the series termwise. 
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Our next result shows that the correlation function © has fractal index a = 2 — 2 q, so that the corresponding 
Gaussian particles have surfaces with Hausdorff dimension 2 + q, as illustrated in Figure [T) 


Theorem 2. If 0 < q < 1, the correlation function © satisfies 


lim 

640 


C(0) - <7(0) 

9 2 - 2 i 


= b 




( 10 ) 


where 

c\ 2,Q P OO /*7T 

b q = —— J xl q J \X~ q — ( x2 + 1 — 2x cos (jf) q “J d cp dx 

Tr 2 q +1 r(i-±q) 2 r(q) 
c 2 (i- q ) 2 r(ig) 2 r(i- 9 )' 

In particular, the correlation function has fractal index a = 2 — 2q. 

We defer the proof of this result to the Appendix. The power kernel © has a negative exponent and thus is 
unbounded, which may lead to unbounded particle realisations. While values of q < 0 are feasible, they are of less 
interest, as the associated correlation functions have fractal index a = 2, thereby generating smooth particles only. 

5. Examples 

Here, we demonstrate the flexibility of the Gaussian particle framework in simulation examples. First, we 
introduce a simulation algorithm. Then we simulate celestial bodies whose surface properties resemble those of the 
Earth, the Moon, Mars, and Venus, as reported in the planetary physics literature. Furthermore, we study and 
simulate the planar particles that arise from the two-dimensional version of the particle model. 

5.1. Simulation algorithm 

To sample from the Gaussian particle model we utilise the property that the underlying measure is indepen¬ 
dently scattered. Specifically, for every sequence {A n ) of disjoint Borel subsets of § 2 , the random variables L(A n ), 
n = 1,2,... are independent and L(UA n ) = Y^L(A n ) almost surely. Let (A n )% =1 denote an equal area partition 
of 8 2 , so that A (A n ) = 47 r/V for n = 1,..., N. The random field X in ([!]) can then be decomposed into a sum of 
integrals over the disjoint sets A n , in that 

JV 

X(u) = J2 / k(v, u)L(dv), uG§ 2 . 

n—i J A n 

For n = 1,..., N fix any point v n € A n . We can then approximate the random field X by setting 

N 

x(u ) = y; k(v n , u ) L(A n ), «e§ 2 . 

71—1 


( 11 ) 

( 12 ) 
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Table 1: Mean radius ro, difference d+ between maximal and mean radius, and difference d_ between minimal and mean 
radius, for Venus, Dry Earth, the Moon, and Mars, in kilometres. 


Body 

Venus 

Dry Earth 

Moon 

Mars 


6051.8 

6367.2 

1737.1 

3389.5 

d+ 

11.0 

8.8 

5.5 

21.2 

d- 

-3.0 

-11.0 

-12.0 

-8.2 


Let us denote the multivariate normal joint distribution of L(Ai),..., L(An) by Fy. To simulate a realisation y of 
the particle Y c , we use the following algorithm. 

Algorithm 1. 

1. Set M = M X M 2 , where Mi and M 2 are positive integers, and construct a grid on § 2 . Using 

spherical coordinates, let u m = and put 0 m = m/M\ and 4> m = 2nj/M 2 , where m = iM 2 + j for 

i = 0,l,..., M\ — 1 and j = 1 , . .., M 2 . 

2. Apply the method of Leopardi (2006) to construct an equal area partition Ai,..., An of S 2 . 

3. For n = 1 ,... ,N, let v n have spherical coordinates equal to the mid range of the latitudes and longitudes 
within A n , respectively. 

4. For n = 1 ,..., N, generate independent random variables L n from Fn- 

5. For m=l,...,M, set x(u m ) = max(c, J2n=i H v n, u m ) L n ). 

6. Set y to be the triangulation of {( u m , x(u m )) : m = 1,..., M}. 


The equal area partitioning algorithm of [Leopardi (.2006 ) is a recursive zonal partitioning algorithm. That is 


after appropriate polar cap areas have been removed, the sphere is divided into longitudinal zones, each of which is 
subsequently divided by latitude. By construction, the equal areal partition cells are continuity sets with respect 
to the intensity of the Gaussian measure, and in all our examples the kernel k is continuous almost everywhere. 


This simulation procedure has been implemented in R (R Development Core Team 


2009), and co de is available 


from the authors upon request. It can be considered an anal ogue of the moving average method (lOliver 


Cressie and I’avlicova. 


2002 


Hansen and Thorarinsdottir 


1995 


20131) for simulating Gaussian random fields on Euclidean 


spaces. In principle, M and N can take any positive integer values. However, the usual trade-off applies, in that the 
quality of the realisations increases with M\, M 2 , and N, at the expense of prolonged run times. For the realisations 
in Figures [l][5l we used Mi = 200 and M 2 = 400, or M = 8 x 10 4 , and N = 10 5 . 

5.2. Celestial bodies 

The geophysical literature has sought to characterise the surface roug hness of the Earth a nd other celestia l bodie s 


in th e solar system via the Hausdorff dimension of their topography ( Mandelbrot 


with 


1982; 


Kucinskas et al 


1992), 


Turcotte (1987) arguing that the dimension is universal and equals about 2.5. Here, we provide simulated 


version of the planets Earth, Venus, and Mars, and of the Moon, under the Gaussian particle model m, with k 
being the power kernel ([5]). We set q = which gives the desired fractal dimension for a Gaussian particle surface, 
and choose the parameters p = rg jci and cr 2 = (d + — d_ ) 2 /c 2 of the Gaussian measure such that they correspond to 
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Figure 2: Simulations of Venus, the Earth, the Moon, and Mars in true relative size. 




Figure 3: Radial function along the equator for the simulated bodies in Figure [2] in kilometres. Clockwise from upper left: 
Venus, the Earth with ocean level indicated by a dashed horizontal line, Mars, and the Moon. 


Price ( 198811 


Jones and Stofan 


reality . For this we use the information listed in Table [T| which was obtained from 
(120081 ). and online sources. The values concerning the Earth describe ‘Dry Earth’; to simulate ‘Wet Earth’ we make 
a cut-off that corresponds to the Gaussian particle Y c with truncation parameter c = 6371 kilometres. In principle, 
we also need to make a cut-off at c = 0 for ‘Dry Earth’, but this is unnecessary in essentially all realizations. In our 
simulation algorithm, we use M\ = 200, M 2 = 400, and N = 10 6 to obtain the celestial bodies in Figure [2j The 
corresponding radial functions along the equator are shown in Figured 


5.3. Planar particles 

We now reduce the dimension and consider the planar Gaussian random particle 

Y c = {o + ru : 0 < r < max(X(u), c)} C R 2 . 
ues 1 
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Table 2: Analytic form, parameter range, constants and associated fractal index for parametric families of isotropic kernels 
k : [0, 27r) —> R on the circle S 1 . 


Kernel 

von Mises-Fisher 

Uniform 

Power 

Analytic Form 

k{6) = e a cos 8 

k{0) = 1(0 < r) 

m = ay q ~i 

Parameter 

a > 0 

re(0,f] 

ge(-i,0)U(0,i) 

Cl 

2nl 0 (a) 

2 r 

2 tt 9 

1-9 

C2 

2nIo(2a) 

2r 

4 7T ^ 

1-39 + 2(72 

Fractal Index 

2 

1 

1-29 


Table 3: Values of the parameter a for the von Mises-Fisher kernel, the parameter r for the uniform kernel, and the 
parameter q for the power kernel used to generate the planar particles in Figure [4] 


Row 

a 

r 

9 

1 

3 

1.5 

0.05 

2 

30 

1.0 

0.25 

3 

300 

0.5 

0.45 


Here c > 0, o G R 2 is an arbitrary centre, and the radial function X(u) is modelled as 


X(u) = f K (u, u) L{ du), u G S 1 , 

J s 1 

with a suitable kernel function K : S 1 x S 1 —> R and a Gaussian measure L on the Borel subsets of the unit sphere 
S 1 = {x G R 2 : ||x|| = 1}. We further construct planar gamma particles where L is a gamma measure on the 
Borel subsets of the unit sphere, L(A) ~ Gamma (k A(A), t) with shape n > 0 and rate r > 0. Under this model, 
/ix = kc\/t and a\ = kc^/t 2 . The Gamma measure is independently scattered and we can thus apply the same 
simulation method as for the Gaussian particles. 

As previously, we assume that the kernel function K is isotropic, in that I\(v,u) = k(d(v,u)) depends on the 
points v,u G S 1 through their angular or circular distance d(v,u) G [ 0 , 7 r], only. Table [2] lists circular analogues of 
von Mises-Fisher, uniform, and power kernels along with analytic expressions for the integrals 


c„ = f k(d(v, u)) n du = 2 f k(rj) n dr], 

J § 1 Jo 


where n = 1,2, and the fractal index, a , of the a ssociate d correlation function, as defined in equation (|3|). The 


power kernel model has previously been studied bylWood ( 1995L Example 3.3). In analogy to the respective result 
on § 2 , if max„ g gi X (u) > c, the boundary of the Gaussian particle Y c has Hausdorff dimension D = 2 — ^ almost 
surely. 
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Figure 4: Planar particles with mean /ix = 25 and variance = 10. Columns 1 and 2 show particles generated using 
a von Mises-Fisher kernel, columns 3 and 4 particles using a uniform kernel, and columns 5 and 6 particles using a power 
kernel, with parameters varying by row as described in Table [3] The particles in columns 1, 3, and 5 are generated under a 
Gaussian measure, those in columns 2, 4, and 6 under a gamma measure. 


The general form of the associated correlation function is 

rir—9 


■ - 0 - 0) d (j>+ / k(<f>)k(0 + d0 

Jo 


C{9) = —( [ k((j>)k(2n - cj) - 9) 

^2 y J 7 r —0 

+ J k((f>)k(9 — <j>) d(f> + J k(<p)k((j) — 9) 0<9<ir. 

For the von Mises-Fisher kernel with parameter a > 0, the correlation functions admits the closed form 

I 0 2(1 + cosd)^ 


C{9) = 


do (2a) 


0 < 9 < 7T, 


where Iq denotes the modified Bessel function of the first kind and of order 0, and for the uniform kernel with 
cut-off parameter r £ (0, f], we have 

C(9) = M — 1(0 < 2r), 0<9<n. 


For the power kernel with parameter q £ (0, 5), tedious but straightforward computations result in a complex closed 
form expression, and a Taylor expansion about the origin yields the fractal index, a = 1 — 2 q, stated in Table [2] 
Thus, the von Mises-Fisher and uniform kernels result in Gaussian particles with boundaries of Hausdorff 
dimension 1 and |, respectively. Under the power kernel, the Hausdorff dimension of the Gaussian particle surface 
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Figure 5: Gamma particles with mean /j,x = 100 and variance a\ = 10, using a gamma measure in © and the power 
kernel © with q = 0.05 (left), q = 0.25 (middle) and q = 0.5 (right). 


is 7j + q. Simulated planar Gaussian and gamma particles with von Mises-Fisher, uniform, and power kernels 
are shown in Figure © with the parameter values varying by row, as listed in Table [3] The simulation algorithm 
of Section o continues to apply with natural adaptions, such as defining the simulation grid u m = 2mn/M for 
to = 1,..., M, where we use M = 5,000 and N = 10 5 . 


6. Discussion 


We have proposed a flexible framework for modelling and simulating star-shaped Gaussian random particles. 
The particles are represented by their radial function, which is generated by an isotropic kernel smoothing on 
the sphere. From a theoretical perspective, the construction is perfe ctly general, as every continuous isotropic 
correlation function on a sphere admits an isotropic convolution root (iZiegel 12014). The Hausdorff dimension of 
the particle surface depends on the behaviour of the associated correlation function at the origin, as quantified by 
the fractal index. Under power kernels we obtain Gaussian particles with boundaries of any Hausdorff dimension 
between 2 and 3. 

While a non-Gaussian theory remains elusive, we believe that similar results hold for gamma particles where 
L(A) ~ Gamma (k A(A), r) in ©) with shape k > 0 and rate r > 0. For instance, Figures ©and [5] show Gaussian 
and gamma particles under the power kernel, respectively. The surface structure for the different bases resemble 
each other, even though the particles exhibit more pronounced spikes under the gamma basis, see also the planar 
particles in Figure © Simil ar parti cle models may be generated using different types of Levy bases L, such as 


Poisson or inverse Gaussian (Jonsdottir et al 


200811 . 


We have focused on three-dimensional particles, except for brief remarks on planar particles in the preceding 
section. However, the Gaussian particle approach generalises readily, to yield star-shaped random particles in 
for any d > 2. The particles are represe nted by their radial function and associated with an isotropic random 
field on the sphere S d_1 . In this setting, [Estrade and I stas (120101 1 derive a recursion formula that yields closed 
form expressions for the isotropic correlation function on that arises under a uniform kernel. In analogy to 


19.93), we refer to this correlation function as the ‘spherical 


terminology used in the Euclidean case (iGne iting. 
hat’ function with cut-off paramete r r g (0, Any spherical hat function has a linear behaviour at the origin, 


and thus has fractal index a = 1. 


Estrade and Istas ( 2010 ) also show that scale mixtures of the spherical hat 
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fu nction provide correlation funct ions of any desired fr actal index a G (0,1], similarly to the corresponding results 


of 


Hammersley and Wider (1955) and 


Gneiting ( 199911 in the Euclidean case. 


A far-reaching, natural extension of our approach uses non-isotropic kernels to allow for so-called multifractal 
particles, where the rough ness proper tie s and the Hausdorff dimension may vary locally on the particle surface. 


This fits the framework of 


Gagnon et al. 


(I2006I1 . who argue that the topography of Earth is multifractal, and allows 


for multifractal simulations of three-dimensional celestial bodies, as opposed to extant work that applies to the 
topography of ‘Flat Earth’. 

We have not discussed parameter estimation under our modelling approach, leavi ng t h is to future wo rk. In a 
Bayesian setting, inference could be performed similarly to the methods developed by Wolpert and Ickstadt (19981), 
who use a construction akin to the random field model in (JTJ) to represent the intensity measure of a spatial point 
process, and propose a simulated inference framework, where the model parameters, the underlying random field, 
and the poin t process are updated in turn, conditional on the current state of the other variables. Alternatively, 


Ziegell (120131 ) proposes a non-parametric inference framework based on series of Gegenbauer polynomials. 


Acknowledgements 

The authors thank Anders Rpnn-Nielsen, Eva B. Vedel Jensen, Jens Ledet Jensen, Richard Askey, Werner Ehm, 
two anonymous reviewers, and the editor for comments and discussions. This research has been supported by the 
Centre for Stochastic Geometry and Advanced Bioimaging at Aarhus University, which is funded by a grant from the 
Villum Foundation; by the Alfried Krupp von Bohlen und Halbach Foundation; by the German Research Foundation 
(DFG) within the programme “Spatio-/Temporal Graphical Models and Applications in Image Analysis”, grant 
GRK 1653; by Statistics for Innovation, sfi 2 , in Oslo; by the U.S. National Science Foundation, grant DMS-1309808; 
and by a Romberg Guest Professorship at the Heidelberg Graduate School for Mathematical and Computational 
Methods in the Sciences, funded by the German Universities Excellence Initiative grant GSC 220/2. 


References 

Adler, R. J. (2010). The Geometry of Random Fields (SIAM Classics ed.). Philadelphia: SIAM. 

Cressie, N. and M. Pavlicova (2002). Calibrated spatial moving average simulations. Statistical Modelling 2, 
267-279. 

Digital Library of Mathematical Functions (2011). Release 2011-07-01, http://dlmf .nist.gov. 
do Carmo, M. P. (1976). Differential Geometry of Curves and Surfaces. Englewood Cliffs: Prentice-Hall. 

Estrade, A. and J. Istas (2010). Ball throwing on spheres. Bernoulli 16, 953-970. 

Falconer, K. (1990). Fractal Geometry: Mathematical Foundations and Applications. Chichester: John Wiley and 
Sons. 

Fisher, N. I., T. Lewis, and B. J. J. Embleton (1987). Statistical Analysis of Spherical Data. Cambridge: Cambridge 
University Press. 

Gagnon, J.-S., S. L. Lovejoy, and D. Schertzer (2006). Multifractal earth topography. Nonlinear Processes in 
Geophysics 13, 541-570. 

















Gaussian random particles 


15 


Gangolli, R. (1967). Positive definite kernels on homogeneous spaces and certain stochastic processes related to 
Levy’s Brownian motion of several parameters. Annales de VInstitut Henri Poincare section B 3, 121-226. 

Gneiting, T. (1999). Radial positive definite functions generated by Euclid’s hat. Journal of Multivariate Analysis 69, 
88-119. 

Grenander, U. and M. I. Miller (1994). Representations of knowledge in complex systems. Journal of the Royal 
Statistical Society Series B 56, 549-603. 

Hall, P. and R. Roy (1994). On the relationship between fractal dimension and fractal index for stationary stochastic 
processes. Annals of Applied Probability 4, 241-253. 

Hammersley, J. M. and J. A. Nelder (1955). Sampling from an isotropic Gaussian process. Proceedings of the 
Cambridge Philosophical Society 51, 652-662. 

Hansen, L. V. and T. L. Thorarinsdottir (2013). A note on moving average models for Gaussian random Helds. 
Statistics and Probability Letters 83, 850-855. 

Hausdorff, F. (1919). Dimension und auBeres MaB. Mathematische Annalen 79, 157-179. 

Hobolth, A. (2003). The spherical deformation model. Biostatistics 4, 583-595. 

Hobolth, A., J. T. Kent, and I. L. Dryden (2002). On the relation between edge and vertex modelling in shape 
analysis. Scandinavian Journal of Statistics 29, 355-374. 

Hobolth, A., J. Pedersen, and E. B. V. Jensen (2003). A continuous parametric shape model. Annals of the Institute 
of Statistical Mathematics 55, 227-242. 

Jones, R. H. (1963). Stochastic processes on a sphere. Annals of Mathematical Statistics 34, 213-218. 

Jones, T. and E. Stofan (2008). Planetology: Unlocking the Secrets of the Solar System. Washington, D.C.: National 
Geographic Society. 

Jonsdottir, K. Y., J. Sclnniegel, and E. B. V. Jensen (2008). Levy based growth models. Bernoulli 14, 62-90. 

Kent, J. T., I. L. Dryden, and C. R. Anderson (2000). Using circulant symmetry to model featureless objects. 
Biomet.rika 87, 527- 544. 

Kucinskas, A. B., D. L. Turcotte, J. Huang, and P. G. Ford (1992). Fractal analysis of Venus topography in Tinatin 
Planitia and Ovda Regio. Journal of Geophysical Research 97, 13635- 13641. 

Leopardi, P. (2006). A partition of the unit sphere into regions of equal area and small diameter. Electronic 
Transactions on Numerical Analysis 25, 309-327. 

Mandelbrot, B. B. (1982). The Fractal Geometry of Nature. New York: W. H. Freeman and Company. 

Miller, M. I., S. Joslii, D. R. Maffitt, J. G. Mcnally, and U. Grenander (1994). Membranes, mitochondria and 
amoebae: Shape models. Journal of Applied Statistics 21, 141- 163. 

Muinonen, K., T. Nousiainen, P. Fast, K. Lumme, and J. I. Peltoniemi (1996). Light scattering by Gaussian random 
particles: Ray optics approximation. Journal of Quantitative Spectroscopy and Radiative Transfer 55, 577-601. 

Munoz, O., H. Volten, J. W. Hovenier, T. Nousiainen, K. Muinonen, D. Guirado, F. Moreno, and L. B. F. M. 
Waters (2007). Scattering matrix of large Saharan dust particles: Experiments and computations. Journal of 
Geophysical Research 112, D13215. 

Oliver, D. S. (1995). Moving averages for Gaussian simulation in two and three dimensions. Mathematical 
Geology 27, 939-960. 

Orford, J. D. and W. B. Whalley (1983). The use of the fractal dimension to quantify the morphology of irregular¬ 
shaped particles. Sedimentology 30, 655-668. 

Price, F. (1988). The Moon Observer’s Handbook. Cambridge: Cambridge University Press. 



16 


L. V. Hansen et al. 


R Development Core Team (2009). R: A Language and Environment for Statistical Computing. Vienna, Austria: 
R Foundation for Statistical Computing. 

Stoyan, D. and H. Stoyan (1994). Fractals, Random Shapes and Point Fields. Chichester: John Wiley & Sons. 

Tovcliigrecliko, A. and I. A. Vakser (2001). How common is the funnel-like energy landscape in protein-protein 
interactions? Protein Science 10, 1572-1583. 

Turcotte, D. L. (1987). A fractal interpretation of topography and geoid spectra on the Earth, Moon, Venus, and 
Mars. Journal of Geophysical Research 92, E597-E601. 

Wicksell, S. D. (1925). The corpuscle problem. A mathematical study of a biometric problem. Biometrika 17, 
84-89. 

Wolpert, R. L. and K. Ickstadt (1998). Poisson/gamma random field models for spatial statistics. Biometrika 85, 
251-267. 

Wood, A. T. A. (1995). When is a truncated covariance function on the line a covariance function on the circle? 
Statistics & Probability Letters 24, 157-164. 

Xue, Y. and Y. Xiao (2011). Fractal and smoothness properties of space-time Gaussian models. Frontiers of 
Mathematics in China 6, 1217- 1248. 

Ziegel, J. (2013). Stereological modelling of random particles. Communications in Statistics. Theory and Methods 42, 
1428-1442. 

Ziegel, J. (2014). Convolution roots and differentiability of isotropic positive definite functions on spheres. 
Proceedings of the American Mathematical Society 142, 2053-2077. 


Appendix A. Proof of Theorem [2] 

We proceed in two parts, demonstrating first the asymptotic expansion m with the constant b q in CD: and 
then establishing the equality of the expressions in CD and CD: which confirms that b q is strictly positive. The 
claim about the fractal index then is immediate from Theorem [lj 

In what follows, if A(-) and B(-) are nonnegative functions on a common domain, we write 

A<B 


if there is a constant C > 0 such that A < CB and C is independent of any parameters or arguments appearing in 
A and B when the latter are allowed to vary in their specified domains. 


A.l. Asymptotic expansion (HOI) 

Recall that 


C{6) = — I (7t 9 A q — l)sinA f (n q a(0, A, 4 >) q — l) d0 dA 
c 2 Jo Ja(\) 


where 


t(0, A, (f >) = sin 0 sin A cos (j) + cos 6 cos A, a{6, A, cf) = arccosf(0, A, cf), 


and 

A(A) = {<p € [0,7r] : 0 < a{9, A, <j>) < 7r}. 
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Therefore, 


^(<7(0) -C{9)) = [ ( 7 r 9 A -9 — l) sin A < f ( 7 t 9 A ” 9 - l) d0 - [ (7r q a(9, A, 0 )" 9 - l) # l dA. 

2 J o [Jo Ja{\) J 

Since A(A) = [0, 7 r] for A G (0, tt — 9] and A{\) C [0, 7 r] for A G ( 7 r —0, 7 t), we decompose the integral on the right-hand 
side as Pi g (0) + P 2 q(d), where P\ q (9) and P2q(9) correspond to the integral with respect to A over (0, 7 r — 9) and 
( 7 r — 9, 7r), respectively. The first mean value theorem for integration implies that there exists a £ G (tt — 6, ir) such 
that 

P2q(0) = 9 [TT q t~ q — l) sint < f [TT q t~ q — l) d</> — [ [i: q a(9 1 t, f>)~ q — l) d<j) >. 

[Jo J-A(t) J 

Hence, -P 2 g (0) decays at least as fast as 0(9 2 ) as 9 { 0. 

As regards the first term, substituting A = 9x yields 


P lq (9) = 9 2 ~ 2q n 2q 


r '( 7r e P B sin (9x) , _ 


(x~ q -ir- q 9 q ) / [x~ q — a(0,9x, f)~ q 9 q ') dcj> dx. 


In order to prove the asymptotic behaviour m it now suffices to show that 

/•OO 

where 

r(r-8)/8 • n (Q x \ 

1(9) = — { —P-(x- q -n- q 9 q )f(9,x)dx 


(13) 


for 9 > 0, with 


f(9,x) — f (x 9 — a(9,9x, </>) 9 0 9 ) d(j> 

Jo 


for x > 0 and 9 > 0. As we aim to find the limit limgj,o /(0), we may assume that 9 G (0, 9q) for some 0 < 9o <C 1, 
and that A G [0, 7 r — 9\. 

Lemma 1 . We have t(9, A, <j>) < cos (9 — A) and 


i(9,\,(j>)- q < \9-X\~ q . 


Proof. We write 

t(9, A, (f) = sin 9 sin A cos <j> + cos 9 cos A = cos(0 — A) + sin 0 sin A(cos (j> — 1). 

Since cos <j> — 1 G [—2, 0] and the inverse cosine function is monotonically decreasing, the claims follow. 
Lemma 2. Define f(0,x) = limg^o f(9,x). Then the limit exists and equals 
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for x {0,1}. 

Proof. For x fL {0,1} fixed, the integrand in the definition of f(9,x) is bounded in <j>. The claim follows 
immediately from the limit 

lim= .(9.fa,0) = jco /2 

<40 9 

along with Lebesgue’s dominated convergence theorem. Indeed, noting that 


arccos (t) arccos (l — y 2 ) y 
6 y 8 


y = (l-t)l/2 


for t G (0,1), we find that 

a(6,6x, 4>) d , 2 , 

inn--- = arccos (1 — y ) 

010 8 d y v y 1 


(1 — cos 8 cos 8x sin 9 sin 8x 


y=o 


lim , „ 9 

0io \ 8 2 


8 2 


1/2 


cos </> ) = ( x 2 — 2x cos <^> + l) 


1/2 


Lemma 3. We have 

for x G [0, (t r - 8)/8\. 

Proof. We hnd from Lemma |T] that 

and the claim follows. 


\f{9,x)\ <ir(x q + \x - 1| 9 ), 


a(9,9x,(f>) q 8 q < |a? — 11 9 , 


For later purposes we need to find the Taylor expansion of a(0, A, </>) q around 0 = 0, 


a(0, A, 4,)~ q = a( 0, A, cf)~ q + —a(0, A, 4>)~ q + R{G, A, <f>), 


where R denotes the error term. 

Lemma 4. We have 

y , (f)~ q = y~ q + 9 + R{&, V-, 

where the error term satisfies 


\R(d,vA)\<9 2 


+ 


\y-9\ q+2 \y — 9\ q+l sin(y — 8) 


for y G [20,7T - 0], 

Proof. From f(0, y, <f>) = cosy, the zeroth-order term is immediately seen to be y~ q . For the first-order term, we 
compute 

-^-a(8,y, fi) = - 1 (cos0siny cos</> — sin0cosy), 

d0 VI - t 2 
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SO 


Hence, 


d0 


a{9,y,(/)) 


0=0 


= — COS ( 




0=0 a(0, y, cf>)i +1 d9 


a(9,y,(j)) 


q COS <f) 


0=0 y 


g+1 


To estimate R, we present it in Lagrange form, 


1_^ 

2~d9 2 ' 


R{0,y,4>) = 7;-TE2a,(9,y,<l>) q 8 2 , 


0=0 


for some 8 £ [0,0]. Then, we have 


d0 2 


i(0,y,<l>) q = 


dd \ a(0, y, </>) 9+1 d0 


a(0,y,4>) 


q(q + l) ( d 


a(0, y, <(>) 9+2 \ d0 


—a(0,j/,(/>) - 


a(0, y, (j>) qJrl d0 2 


a(0,y,</>) 


= A^iy^) + B{9,y,(j)). 


Let us now estimate H(0, y, </>) from above for y £ [20,7r — 0]. From 


< 


we get 


\/l — t 2 y/l — cos 2 (y — 0) sin(y — 0) 
sin y -f sin 0 




< 


sin(y — 0) 

By monotonicity, we have sin0 < sin(y — 0) and siny < 2 sin(y — 0). Hence, in view of Lemma [U we find that 


\M9,y,<!>)\ ;$ \y-o\ 


-9-2 


for y G [20,7r — 0]. 

To estimate H(0, y, <^) in the same range, we compute 


—a(0, y, 0) = (7(0, y, 0) + D(0, y, 0), 


where 


and 


(7(0, y,<j>) = —— (- ) (cos0sinycos 4> — sin0cosy) = —7- ^ . (cos0sinycos <j) — sin0cosy) 2 

d0 V v 1 -t 2 / (l-t 2 )^/ 2 


—-— (cos 0 sin y cos </> — sin 0 cos y) = 1 (sin 0 sin y cos <j) + cos 0 cos y). 


^[T^^ 2 


D(0,y, <j>) = 


VT^f- d 0 
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We have 


for y G [26, n — 9\. Hence, 


1*1 


< |cos (y - 6 )I < 1 


(1 — t 2 ) 3 / 2 sin 3 (y — 9) sin 3 (y — 9) 


i/ (siny + sin6>) 2 1 

|C(6>,y,(/))| < --- < 


sin 3 (y — 6) "sin (y - 9)' 


Similarly, \D(8,y,<j>)\ < l/(sin (y — 9)). Hence, 

\B(9,y,(ft)\ < 

for y G [29, ix — 9\. 

Finally, for 9 G [0, 9\, 


\y - 9\ q+1 sin(y - 9) 




< 1 


and similarly, 


B(9,y,<j>) 


< 


(y-9)*+ 2 " |y-0r +2 ’ 

i_<_ 1 


(y — 9) q+1 sin (y — 9) ~ \y — 0| 9+1 sin(y — 9) 
Combining the estimates for A and B, the proof of the lemma is complete. 


In what follows, we need the classical estimate 


sin 9x 

0 < —-— < x, 


(14) 


for x > 0 and 9 > 0. 


Lemma 5. We have 


forx G [2, (n — 9)/9]. 


sin 9x 
9 


l/(M)| < (x- 1) 


-i - 9 


Proof. From Lemma 0] and cos (ft dtp = 0, we find 


f(9, x) = 9 q+2 f R{9,9x, (ft) d<ft. 
Jo 


Using Lemma 0] and (TUH) . we get 


5^ 1/(9,x)| < 2^(x - 1)-’-. + - l)-« < (x - 

6 0 smd[x — 1) 


Lemma 6. We have that x 1 q f( 0,x) is Lebesgue integrable and 

pOO 

lhn/(0) = J x 1 ~ <? /(0, x) dx. 
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Proof. From (I14[) . Lemma [3l and Lemma [5[ we find that 




< 


X 1 2q + \x — 1| q , 

(x - l)- 1 - 2 ? + (a; - l)” 1 " 9 , 


x G [0, 2] 
x G [2, oo) 


uniformly in 9 G (O,0 o ], where 0 < 6* 0 <§G 1- Since the latter function is Lebesgue integrable, the claims follow from 
Lebesgue’s dominated convergence theorem along with Lemma [5] 

This completes the proof of (IT51) and therefore of the asymptotic relationship (flUl) with the constant b q in (fill) . 


A.2. Equality of the expressions in (fTTl) and (fl2l) 

We now show that the constant b q is strictly positive. Specifically, we demonstrate the equality of the expressions 
in m and (fl2l) for q G (0,1). Toward this end, we first prove that 


7T 2 « +1 r(l - \q) 


bq c 2 T(ig) J Q 


j ^ t q 1 (1 -e * ifi(l - \q\ l;t)) 


(15) 


where with (x)n = 1 a nd (x) n = x(x + 1) • • • (x + n — 1) for n = 1 , 2 ,..., the classical confluent hypergeometric 


function (Digita l Lib rary of Mathematical Functions 


2011 


iFi(a;b; t) = ^ 


k=0 


Chapter 13) can be written as 

(a) k t k 

(b) k kl ' 


We establish this representation as follows. With a keen eye on the inner integral in dm we note that for x > 0 
and </> G (0,7r), 


x~ q = ( x 2 )~ q/2 = 


r (k) 


3 -tx 2 t q/2 fjj 


2 D J 0 


and 


(1 + x 2 — 2xcos <j>) 


-9/2 _ 


i r 
ih) Jo 


t ’ 


a —1( l-\-x 2 —2x cos 4>) -f-q/2 


d t 

T* 


r(|«) jo 

Substituting these formulae into m, and interchanging the order of the integration with respect to (f> and t , we 
obtain 


, 2n 2q 

^9 r /l \ / 

C2T( 1 q) Jo 


j x 1 " 9 J t q / 2 e~ tx2 j (i - e - t(1 ” 2a:cos?i) ) d^ y dx. 


(16) 


By well-known formulae, 


[ e tcos0 d0 = irloit) = tt 0 Fi(1; \t 2 ), 

Jo 


where Iq denotes the mo dified Bessel function of th e first kind_of order 0, and oFi is a special case of the generalised 

formulae 10.32.1 and 10.39.9). Therefore, 


hypergeometric series (Digital Library of Mathematical Functions 


2011 


J (l-e- t(1 - 2xcos ^) d0 = 7r (l — e _t qFi(1; t 2 x 2 )). 
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Substituting this result into (1161) . and interchanging the order of the integration, we obtain 


2 tt 2 « +1 


bq c 2 r(ig) J 0 u J 0 


t ,/2 


x 1 ~ q e~ tx 


(l — e 4 0 -Fi(l; t 2 x 2 )) dx 


d t 


(17) 


With the substitution x = u 1 / 2 , we get 


roo -| 

/ x 1 ~ q e~ tx2 dx = - 12 g - 1 r(l — ig). 
Jo 2 


2011 


We apply next a well-known formula for the Laplace transforms of generalised hypergeometric series ((Digital Li brar y of Mathemg 
formula 16.5.3) to obtain, with the substitution x = u 1 / 2 , 


r OO 1 ^ 

/ a: 1-9 e _4x2 oi 7 i(l;t 2 £ 2 ) da; = -t2 9_1 F(l- ±g)iFi(l- ig;l;i). 
Jo Z ~ 


Consequently, by m, we have established the representation (fl5l) . 


Finally, we apply the Kumrner formula for the \F\ function (Digital Library of Mathematical Functions 
formula 13.2.39) to show that 


2011 


1 — e 4 iFi(1- iqr; 1;t) = -qt 2 F 2 (±g + 1,1; 2,2; -t). 

Thus, by the representation 

n 2q+1 q T(1 — ^q) r°° 


bn = 


2c 2 TUq) j 0 


r°° dt 

/ t q 2^2(|? + 1) 1; 2, 2; — t) —, 
Jo f 


and this integral is a well-known Mellin transform; see 


Digital Library of Mathematical Functions (2011 


16.5.1), where the integral is given in inverse Mellin transform format (Digital Library of Mathematical Functions 


2011 


formula 


Section 1.14(iv)). The proof of the equality of the expressions in (fill) and (fl^l) is now complete. 
















